Draft version May 31, 2013 

Preprint typeset using I^'T^]X style eniulatoapj v. 5/2/11 



ON THE JITTER RADIATION 
S.R. Kelner 

National Research Nuclear University (MEPHI), Kashirskoe shosse 31, 115409 Moscow, Russia; Max-Planck-Institut fiir Kernphysik, 

Saupfercheckweg 1, D-6917 Heidelberg, Germany 



o 

(N 



o 

m 
O 



F.A. Aharonian 

Dublin Institute for Advanced Studies, 31 Fitzwilliam Place, Dublin 2, Ireland; Max-Planck-Institut fiir Kernphysik, Saupfercheckweg 1, 

D-6917 Heidelberg, Germany 

D. Khangulyan 

Institute of Space and Astronautical Science/ JAXA, 
3-1-1 Yoshinodai, Chuo-ku, Sagamihara, Kanagawa 252-5210, Japan 
Draft version May 31, 2013 

ABSTRACT 

In a small scale turbulent medium, when the nonrelativistic Larmor radius i?L = mc^ / eB exceeds 
the correlation length A of the magnetic field, the magnetic bremsstrahlung of charged relativistic 
particles unavoidably proceeds in the so-called jitter radiation regime. The cooling timescale of parent 
particles is identical to the synchrotron cooling time, thus this radiation regime can be produced with 
very high efficiency in different astrophysical sources characterized by high turbulence. The jitter 
radiation has distinct spectral features shifted, compared to synchrotron radiation, towards high 
energies. This makes the jitter mechanism an attractive broad-band gamma-ray production channel 
which in highly magnetized and turbulent environments can compete or even dominate over other 
high energy radiation mechanisms. In this paper we present a novel study on spectral properties 
of the jitter radiation performed within the framework of perturbation theory. The derived general 
expression for the spectral power of radiation is presented in a compact and convenient for numerical 
calculations form. 
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1. INTRODUCTION 

Charged particles moving in electric and magnetic 
fields experience effective energy losses via radiation. Be- 
cause of high conductivity, the electric fields in astrophys- 
ical plasmas are typically screened, thus the radiation is 
dominated by interactions with the magnetic field due to 
the so-called magnetic bresstrahlung. The latter is one of 
the major nonthermal radiation processes in astrophysics 
and operates with high efficiency in a large variety of as- 
trophysical environments. In the case of a regular mag- 
netic field or a chaotic field characterized by large scale 
fiuctuations, we deal with the so-called synchrotron radi- 
ation. This process and its implications in astrophysics 
have been studied in great detail (see e.g. Ginzburg & 
^yrovatskii 196£; Rybicki & Lightman 1979). In highly 
turbulent environments, namely, when the nonrelativis- 
tic Larmor radius i?L = mc^/eB does not exceed the 
char actoriotic ccalo of turbulence A, the radiation pro 



exp [—{uj/uJc)], where 



3j'^eB/{2mc) is the 



ceeda in significantly different regimp, which in the as- 



trophysica l literature is referred as diff usive sync hrotron 



radiation (Toptygin fc Fleishma: 



or as jitter radiation ( Medvedc 



!^|1987|, hereafter [TFSTI ) 
^ pOOq , hereafter pro|) . 



Hereafter we will use the term "jitter" . 

The spectral features of jitter radiation substantially 
differ from the synchrotron radiation. While the power 
of the synchrotron radiation of a monoenergetic parti- 
cle Ptj is described with a good accuracy as uiPu, oc 
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characteristic synchrotron frequency, in the case of jit- 
ter radiation the peak is shifted towards higher frequen- 
cies by a factor of a = i?L/A. Unless the distribution of 
emitting particles is strictly monodirectional, the power 
spectrum of jitter radiation below the maximum is fiat, 
i.e. the Spectral Energy Distribution (SED)F] ujP^ cx w , 
while beyond the cutoff energy it has a power-law be- 
havior, uPi^ cx u!^~°' , where a is th e power-law in- 
dex of the turbulence spectrum ( TF87 ). Thus, instead 
of the typical exponential cutoff in synchrotron spec- 
trum, the jitter mechanism yields a power-law spectrum 
which can be extended up to the frequency of a^LOc ■ 
This makes the jitter radiation of electrons an excellent 
high energy gamma-ray production process in contrast 
to the synchrotron radiation which even in the case of 
extreme accelerators operating at the maxi nium possi- 
ble rate allo wed by classical electrodynamics ( Aharonian 



et al. 2002) is limited by the maximum possible energy 
eo = huja = 9/4a~iTOc2 - 150 MeV. 

However, so far this remarkable feature of jitter radia- 
tion practically has not been explored for interpretation 
of high energy gam ma-ray phenomena (see h owever, the 
recent papers by Teraki & Takahara 2013). Instead, 
more emphasis has been placed on the energy int erval 
below the cutoff. In particular, it is augured in MOO 



^ The so-called Spectral Energy Distribution or SED is deter- 
mined as Lj^dN/dw or e^dN/de, where dN/dw (dN/de) is the dif- 
ferential distribution over frequencies (energies). Obviously here 
the SED is uPu ; note that in astronomical literature for SED is 
often used the denotation i/Fu . 



that the jitter radiation below the cutoff can result in 
harder spectra than the synchrotron radiation, namely 
ujPuj oc 10^ . However, the claimed energy dependence is 



close ly related to the assumed geometry of the magnetic 
field. Namely, it can be achieved if the magnetic field 
has only one non-zero component which can be realized 
only for a rather unrealistic configuration of the turbu- 
lent field (see Sect. for details). 

In the seminal paper on jitter radiation by [TFS? it 
has been realized that the spectral maximum of the jit- 
ter emission is located at higher frequencies than in the 
synchrotron regime, and that the high energy part of 
the jitter spectrum could be described by a power-law. 
Thus, even for the case of monoenergetic particle dis- 
tribution one may expect a broken power-law spectrum. 
This should lead to the modification of the standard re- 
lations between spectral slopes, flux levels and breaks 
found in synchrotron spectra. Possible applications of 
the jitter mechanism also has been discussed, basically 
in the low energy range of cosmic electromagnetic radi- 
ation. In particular, it has been proposed that the jitter 
radiation can be responsible for the radio to optical (X- 
ray) s pectra of some active galaxies and pulsar wind neb- 
ulae (^^FSllFleishman fc Bietenhol^|2007t Mao & Wand 

poo7|) ! — 

The underestimation of the potential of jitter radiation 
for production of high and very high energy gamma-rays 
could be related to the effect of weakening of the diffu- 
sive shock acceleration process in the case of short-length 
scale turbulence. A self-consistent consideration of the 
processes of particle acceleration and emission in the 
framework of the diffusive shock acceleration paradigm 
predicts a shift of the jitter radiation peak towards low 
frequencies as comp a red to the pure synchrotron radia- 
tion ( |Derishev||2007[ |Kirk fc Revlu^ ^OlOl) . However, ff 
the inhomogeneities responsible for particle acceleration 
and emission are different, e.g. when these processes oc- 
cur in spatially separated regions, the spectral maximum 
would be shifted towards higher energies making the jit- 
ter radiation a very effective high energy gamma-ray pro- 
duction mechanism. Therefore the spectral features of 
this radiation in the entire energy range deserve detailed 
qualitative studies. 

To explore the process in a general form, we propose 
a new approach based on the perturbation theory. In 
terms of additional assumptions, the proposed method 
is less demanding compared to previous studies, and al- 
lows a precise control of the applicability conditions for 
the derived solutions, e.g. the range of the high energy 
power-law extension beyond the spectral maximum. 

In this regard we should note that in previous studies 
some principal results have been obtained under specific, 
altho ugh not always obvious assumptions. For example, 
M00| has derived the spectrum of radiation for the case 
of a very specific geometry of the magne tic field fluc- 
tuati ons. In som e others studies (see e.g., Fleishman fc 
Bietc nholz 2007) the jitter radiation spectrum in fact has 
not been strictly derived, but rather predefined through 
its asymptotic behavior. Finallv. some studies address 
the c ase of anisotropic turbulence ( Reynolds fc Medvedev 
however the structure of the used correlation ten- 



poT2| ) 



sor IS not consistent with the fundamental requirement 
oi V ■ B — Q. We discuss these concerns in detail in 
Section ^ 



Finally, we should mention that a significant progress 
recently has been achieved through numerical co mputa- 
tions based o n particle-in-cell techniqu e (see, e.g., Reville 
fc Kirk||20To|;|Teraki fc Takahara||201l|). This method has 



a great potential to deal with quite complex distributions 
of emitting particles. On the other had, the analytical 
approach allows better understanding and interpretation 
of physics behind the obtained results. In this regard, 
two methods are complementary and equally important. 
The paper is organized as following: in Section g the 
basic results on the energy spectra, as well as the ap- 
plicability limits for the derived spectra are presented, 
in Section we consider the case of chaotic magnetic 
field. In Section ^we compare the radiation properties in 
chaotic magnetic field with the conventional synchrotron 
radiation; the latter is is briefly discussed in Section Q. 
The case of anisotropic turbulence (under assumption of 
isotropic distribution of emitting particles) is considered 
in Section pi Finally, we compare our results with previ- 
ous studies in Section M, and summarize the main results 
in Section g. 

2. PERTURBATION THEORY 

The intensity and the energy distribution of radiation 
produced by a particle of a given charge e depends only 
on its trajectory. Let r{t) and v{t) — r{t) be the radius- 
vector and the velocity of the particle at the instant t . 
Then, the energy sp ectrum of radi ation is described by 
equation (14.65) of Jacksorj (1998) which for our pur- 
poses is convenient to present in the form|j 



dSr, 
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4,^2^3 



U{t) dt 

j U{tx)U*{t2)dtxdt2. 



(1) 



Here the integrand 



jj,.. _ nx [{n-f3{t)) ^a{t)] ^^^^^) 



(2) 



depends on the particle velocity v{t) — c(3{t) and the 
acceleration a{t) = cf3{t) , as well as the function $(i) = 
Lu{t — nr{t)/c), where n is the unit vector towards the 
momentum of the radiated photon. The function U*{t) 
is the is complex conjugation of U{t). Note that equa- 
tion (|l|) is precise; it is derived within the framework 
of classical electrodynamics through integration of the 
Maxwell equations in vacuum. 

It is convenient to introduce new variables of inte- 
gration: t = {ti + t2)/2 and t = ti — t2 (note that 
dti dt2 = dtdr). Then equation (nl) results in 

dS„,u, 



dujdQ Att^c^ 



f U{t + T/2)U*{t-T/2)dtdT. (3) 

As it is shown bellow, for a fixed value of t the integrand 
rapidly decreases with the increase of \t\ . Also, the inte- 
grand is characterized by a weak dependence on t . The 

^ d£'ni^ is the energy radiated by a particle into the solid angle 
dQ within the frequency interval dtu . 



integration of the integrand over dr gives the spectral 
power of emission at the moment t : 



Pr^Ut) 
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/ U{t + T/2)U*{t-T/2)dT. 



(4) 



Note that if one considers equation (0) as a classical 
limit of the corresponding quantum relation, then the 
integrand in equation (0) can be interpreted as the emis- 
sion probabilitv raultipl ied to the photon energy fiuj (see 
Berestetskiiet al.||l989| , § 90). 

The radiation detected by an observer is produced by 
ensemble of particles occupying a certain region in space. 
We will consider the case of a chaotic magnetic field, 
assuming that statistically averaged (over time or space) 
magnetic field (B) = . 

Let A be the correlation length of the magnetic field. 
If the distance between two chosen points at ri and r2 
exceeds A, then the corresponding magnetic fields Bi 
and B2 can be treated as statistically independent, thus 
the time-averaged product of these fields {BipB2a) = 

{Blp){B2a)^0. 

To obtain the radiation spectrum, the integrand in 
equation (||) should be averaged over all possible con- 
figurations of the magnetic field. It is convenient to per- 
form this procedure in the framework of perturbation 
theory. The acceleration of particle is proportional to the 
strength of the magnetic field B , a = e{f3 x B)/{mj). 
In the first approximation, all other relevant parameters 
can be treated as in the absence of the magnetic field, 
i.e. f3{t±T/2)= f3{t), r{t±T/2) = r{t)±(3{t)T/2. The 
applicability of the approach is discussed bellow. This 
approximation results in 

„2 
Pr,u{t) 



47r2c3(l-n/3)2 



a+a 



(ria+){na^) 



72(1 - n/3)2 



^Ml~rvP)r^^^ (5) 



where a± = a{t ± t/2) and /3 = ^(t) ); 7 = 1/^1-/32 
is the particle Lorentz factor. The derivation of equa- 
tion (|) was performed using the formula for the standard 
double vector product: ax {bx c) ~ b{ac) — c{ab) and 
taking into account that in the magnetic field the acceler- 
ation and velocity vectors are orthogonal. Equation (|5|) 
represents the first non-vanishing term in the expansion 
of the emission spectrum in powers of the magnetic field. 
Our ultimate aim is to derive the emission spectrum 
integrated over the emission angles of photons and aver- 
aged over the magnetic field fluctuations. It is convenient 
to select the z axis to be parallel to the particle veloc- 
ity /3 , and start with averaging over the azimuthal angle 
(j) in respect to the direction of the particle velocity /3 . 
Then, the scalar product of the vectors n and /3 does 
not depend on the azimuthal angle (p, {nf3) = l3cos6. 
Given that a± _L f3 , one obtains 



{{na-^-){na^)) ~ — (a+a_) sin^i 



(6) 



a(t - r/2) ait) a{t + t/2) 




Fig. 1. — Schematic description of the basic geometry adopted 
for computations. Thin solid line represents a segment of the tra- 
jectory of the charge particle, which in the framework of the per- 
turbation approach can be taken as a straight line. The particle 
acceleration is orthogonal to the velocity (which in its turn is par- 
allel to the trajectory). Also it is assumed that the acceleration is 
statistically independent for distances along the trajectory which 
exceed the magnetic field correlation length A. This corresponds 
to the condition \t\ > A/c. 

and, after averaging of equation (ph over , we have 
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1- 



47r2c3(l - n/3)2 ^ 272(1 - n/3)2 ^ 

00 

X /"(a+a_)e""(i-"^)"dT. (7) 



In equation (Q) the charge velocity, (3 ~ l3{t) , is 
treated as a constant (independent of r). The aver- 
aging over the magnetic field configurations results in 
appearance of a correlation function, {a^a_) , under the 
integral. Note that for /3 = const, the acceleration and 
magnetic field have the same statistical properties. In 
particular, (a(i)) = and {a{t + t/2) a{t - t/2)) = 
{a{t + t/2)) {a{t — t/2)) = if the distance between the 
corresponding points exceeds A (i.e., if cBt > A). This 
feature of (a+a-) is illustrated in figured. 

Since the radiation of ultrarelativistic particles is 
strongly beamed towards the direction of motion (9 ^ 
1/7), we will consider only the region of small angles, 
9^1. This allows significant simplifications of calcula- 
tions which result in 

e2 ~4Q , 4^4^ /. 
Pr^Ut) - -^ W .LJ / («+«-) e'"' dT , (8) 



7r2c3 (1-^72612)4 



where 



2^(1+^^)- 



(9) 



Since the above results are derived within the frame- 
work of the perturbation approach, it is iniportant to 
study the range of applicability of equation (pi). 

The integrand in equation (|8|) rapidly decreases in 
the range of |r| > A/c. Therefore, the obtained ex- 
pression describes correctly the emission power if the 
terms neglected at derivation of equation (ph are small 
for I'''! ^ A/c. In the precise equation (^, the denomi- 
nator contains a term d± = 1 — nf3{t ± t/2) . For small 

values of T, we have d± — {1 — n(3) =F n(3 -J . For ul- 
trarelativistic particles the angle 9 between n and /3 is 
smaU (^ 1/7), thus, given the orthogonality of $ and 



/3 , the last term in d± can be estimated as 



eBX 



Since 



this term was neglected at derivation of equation (K 
must be small compared to the first term in d±, w 



it 
lich 



is estimated as ^ '^l"f^ ■ Thus, the range of apphcabihty 
is determined by the condition 

eB\ , , 
^<1. 10 

Furthermore, the exponential term in equation (0) con- 
tains a function, iA = i($(f + r/2) - $(t - r/2)J^ The 
Taylor expansion of the function A gives A = w(l — 
n/3)r — a;n,^r'^/24. In derivation of equation (^ only 
the first term in this expansion has been kept, therefore 
the applicability can be reduced to the condition of ne- 
glecting the second term. Since the function A is in the 
exponent, the condition is a;n/3 (A/c)'^ ^ 1 . The mod- 
ule of the particle velocity in magnetic field remains con- 
stant, /3^ = const, thus \ ^0^ = {$)'^+f30 = 0. Since 
emitted photons and the particle velocities are nearly 
parallel, n k, (3 , the term (n/3 ) can be estimated as 
(n/3 ) w {(3(3 ) = — (/3 )^ . This gives the second condition 
of applicability of equation (^): uj{f3)'^{X/c)^ < 1. By 

expressing the acceleration through the magnetic field 
strength, the condition of applicability of equation (@) 
can be written in the form 



uj<$: 






(11) 



Note that for a homogeneous magnetic field A = oo , 
therefore the standard synchrotron spectrum cannot be 
derived in the framework of perturbation theory. 

Equations ( [lO|) and ( [ll| ) as conditions of applicability 
of the perturbation approach can be interpreted in the 
following way. If a charged particle travels in a region 
filled with magnetic field a path A which is shorter com- 
pared to the trajectory curvature i?, then the particle 
is deflected by an angle 56 ~ X/R. The first condition 
given by equation ( piO| ) implies that 60 <^ I/7. Con- 
cerning the second condition given by equation (f^), it 
is equivalent to the requirement that the segment of the 
trajectory of length of order ~ A, can be treated as a 
straight line. 

The characteristic frequency of the radiation in this 
regime, ujj , can be estimated from first principles. 
Namely, while the emission is formed during the time 
interval (5trad ^ A/c, it is registered during Jiobs = 
6trad (1 — fifS) ~ Stj-iid/l'^ ■ Th us the characteristic fre - 
quency is estimated as (see e.g. Landau & Lifshitz 197a) 



1 



Ulj 



Stoh 



C7 



(12) 



Note that this frequency ujj is independent on the mag- 
netic field strength B . 

It is interesting to compare the characteristic fre- 
quencies, at which the bulk of radiation is produced, 
in highly turbulent and homogeneous magnetic fields 
corresponding to the jitter and synchrotron radiation 
regimes. The characteristic synchrotron frequency can 
be expressed through the nonrelativistic Larmor radius 
Rl = mc^/eB: 



3cY 



3 A 



(13) 



2Rl 2 Rl ' 

It is convenient to express equations ( p^ ) and (^ also 



through Rl : 



^«' 



U ^ Uj 



^J 



(14) 



When these conditions are satisfied, the ratio ujjIujc. ^ 
Rl/X ^ 1, i.e. the characteristic energy of photons 
emitted by charged particles in highly turbulent mag- 
netic field may significantly exceed, by a factor of Rl/X, 
the characteristic energy of synchrotron photons emit- 
ted by same particles in a regular magnetic field of same 
strength. 

Finally, one should mention another constraint on ap- 
plicability of equation (H) related to the plasma effects. 
The basic equation (0) describes emission in vacuum ne- 
glecting the impact of the surrounding plasma. If the 
radiating particle is located in plasma, the latter in the 
frequency range w ^ Wp can be treated as a medium 
with dielectric permittivity e(w) = 1 — cjp/cj^ , where 



I A^ne^np 



(15) 



is the plasma frequency ( n^ , me and e are the number 
density, mass and charge of electrons, respectively). At 
eiijj) « 1, the term (1 - n/3) « (A + 6i^)/2 in ah above 

derived formulas should be replaced by the one corrected 

2 

for the dielectric permittivity, (1 — ^/en(3) sa (^ H — | + 

6'^)/2 . Thus, the influence of the medium can be ignored 
for sufficiently high frequencies, w ^ ^p1 ■ Note that 
the particle Lorentz factor 7 and the plasma frequency 
enter the equation in the form of the combination 1/7^-1- 
{(jjp/uj)'^ ■ Therefore, the influence of the medium can be 
taken into account if in all above equations we replace 
the particle Lorentz factor 7 to 



7*(w) = 



7 



(16) 



^1 + (wp7/w)2 ■ 

(see e.g. Ginzburg fc Syrovatski: 1969 ; Fleishman 200(j ). 
However, as long as we are interested in high frequency 
range uj ^ Wj,7 , for the sake of simplicity we will ignore 
(unless otherwise is stated) the difference between 7 and 



r 



3. DEALING WITH CHAOTIC MAGNETIC FIELD 
The integrand of equation (tel) contains the term 



{a+a-) 



jYi'2^2 






{Spa 



{f3xB+){f3xB^) 



VpVcr) -B+p-B-cr , 



(17) 



which should be averaged over different configurations of 
the magnetic field. Here v = (3/\(3\ is unit velocity vec- 
tor. The magnetic field values B+ and B- corresponds 
to the points where the charged particle is located at time 
instants (i±T/2),i.e. B± ^ B{r{t)±f3{t)T/2,t±T /2). 
The statistical averaging of this expression will result in 
appearance of the correlation function: 



Kp„ = {Bp{riM)Ba{r2M)) ■ 



(18) 



which is a second order tensor. Here, under the statistical 
averaging we suppose a general standard procedure; it 
could be a space-time homogenization or an integration 
ove r an ensemble of field co nfigurations (see e.g. § 118 
of p^jandau fc Lifshitz 1980). Here we assume that the 
field is statistically homogeneous and stationary. This 
implies that the correlation function depends only on the 
difference of the coordinates (ri— r2) and times (^1—^2), 
i.e. K„„ = Knrriri — r2,ti — 12). In this case, (B^) = 



^pa 



Kpp{0) — const. 

It is convenient to present the correlation function Kp^ 
in the form of a Fourier integral: 



KpAr,t)^jKpAq,>f)e^^'i^~-'^^ 



(Pq dK 



(27r)3 2^ ■ ^^^^ 

Since the magnetic field is divergence free (VJ3 = 0), 
Kpcr should satisfy the following conditions 



dKpa-fdxp = , dKpa/dxa = , 



(20) 



which for the Fourier transform Kp^ take on form [the 
transversality condition) : 

Kp^qp = , Kp,q, = . (21) 

While in Section [^ we will briefiy discuss different ten- 
sor structures of the correlation function, here we con- 
sider the case of isotropic turbulence. This results in 
the followi ng fo rm of the correlation function (see e.g. 
Fleishman] ^OOSD ): 



Kpa{q,x) 
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2 i'^- - — 



^{\q\,>c){B^ 



(22) 



Here the constant factor (-B^) is introduced which allows 
^ to meet the normalization condition 



, , d'^q dx 1 
^'^' ^^ 72^ 2^^ " 4^ 
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1. 



The tensor structure given by equation (p2), obviously 
satisfies the transversality condition of equation ([2l|). 

The averaged values of {a+a-) can be expressed 
through the correlation function as 

{a+a-) = {5pa - VpVa) Kpa{c(3T, t) . (24) 

Here we took into account that r_|_ — r_ = c/3t, 
and replaced in the numerator 0^ to 1. From equa- 
tions (|l|) and (H) we find 



((a+a_))e'"^dr = 



2j^2^2 



{B') 



1 



{uqf 



*(g,>^)e*("'5'^-^+") 



^ d'^q dx 



dr. (25) 



q- J (27r)3 27r 

After substitution of /3 by the velocity unit vector 
i^, and integration over dr resulting in a (5 -function 
27r 5(cqf — >f -f- a)), the integral over dx can be com- 
puted analytically: 



' dT 



2m^7^ 



IB' 



1 



(yq) 



^{q,uj -f cqv) 



d^q 

(2^ 



(26) 



Thus, we arrive at the following expression for the energy 
and angular distribution of radiation per unit time 



Pr.u,{t) 
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2, 7^(1 + 7^0 



1 



(^g) 



\I'(g,(I' + cqi> 



(1+726*2)4 
d^q 



(27) 



r J ■ (27r)3 

where w is determined by equation (^. 

Let's consider now the case of steady turbulence, i.e., 
when the correlation function given by equation (|l8| ) 
is time-independent. Then the Fourier image of the 
correlation function contains a ^-function, ^fq. k] = 
^(q) 27r(5(>f), and the normalization condition ( |2^ ) be- 
comes 



*(?) 



d^q 

(2^ 



00 

±^j^,{q)q^dq = l. 



(28) 



Note that the function \E' determines the spectrum of the 
energy density of the stochastic magnetic field, since 



{B' 



{B^ 
16 7r3 



■^iq)q^dq. 



(29) 



In the case of the stationary turbulence, a 5 -functional 
factor. 27rJ(iI; -|- cuq) , appears in the integrand of equa- 
tion (p?!). This makes the integration over dilg (note 
that d'^q = q^dqdVtq ) rather trivial: 
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1 + 



c'^q'^ 



'^{q)qdq. 



(30) 



i/c 



Now we can conduct analytical integration over the 
emitting angles of radiation. The major contribution to 
the integral comes from range of small angles 6 < I/7; 
the contribution from large angles, 6 ^ I/7, is negli- 
gibly small. Thus, applying the standard approach for 
calculations of radiation of ultra-relativistic particles, one 
can adopt dfl — 2t:9 dO , and perform integration over 
from zero to infinity. It is also convenient to introduce a 
new integration variable C = 7^6*^ and change the order 
of integration over C, and q . After performing a trivial 
integration over ^ , we arrive at 



P^t) = 



e^B') 



uiO^{q)qdq, (31) 



u/{2c^^) 



where ^ = 2qcj''/Lu, and 



,^, ^ 3 4 31ne 



(32) 



In the range of integration over dq, the variable ^ 
alters from 1 to 00; while the function u{^) increases 



6 



monotonically from m(1) = to u(oo) = 1. Adopt- 
ing ^ as the integration variable, equation ( |3l| ) can be 
presented in the form convenient for numerical compu- 
tations: 



PUt) 



(B' 



GTr^m^c^ 



uiO 



2C72 



* 



2C72 



f , (33) 



Equation ( |3l| ) is an integral function which depends on 
the turbulence spectrum. However, it obeys some general 
properties not affected by the turbulence. In particular, 
from equation (31) follows that independent of ^(g) , the 
radiation spectrum Pu){t) is a monotonically decreasing 
function of w . This feature becomes obvious after the 
differentiation over w : 



dPu, 



j/(2c72 



du! d^ 



^{q)qdq. 



(34) 



Here it is taken into account that the contribution to the 
derivative from the lower integration limit is null (given 
that u{l) = 0). Since du{^)/d^ > 0, ^(q) > and 
d^/duj < , the integrand is negative, and the integration 
results in dPuj/du! < 0. Thus, this function achieves its 
maximum value at uj = 0, i.e.. 



PUt) < Po 



e\B^) 



■^{q)qdq. 



(35) 



Of course, this estimate is meaningful only if the integral 
in the right side of equation converges. 

Note that the photon e nerg y and the particle Lorentz 
factor enter to equation (|3l[) only in the form of ratio 
w/7^ . Thus, the spectrum P^ is, in fact, a function 
of one argument uj /ujj (for intermediate calculations we 
drop, just for simplicity, the argument t): 



P..P(^]^P{^ 



(36) 



Here P{uj/ljJj) is a monotonically decreasing function. 
In case of absence of other characteristic frequencies in 
the physical setup except for ujj , in the range to <C ujj 
function P is nearly constant, P k, Pq. However, at 
very small frequencies the surrounding plasma may sig- 
nificantly change the behavior of P^j ■ The substitution of 
7 by 7, (in accordance with equation (Uq)) leads to equa- 
tion (|31) in which S, should be replaced by ^, = 2qc'^1/Lu . 
Thederivative 



—- = -2qc-f —^ 
ou) (w^ ■ 



I'^l 



■ 7^^2)2 



(37) 



has a positive sign at w < ^Up , and becomes negative 
when u > ^ujp . Therefore, independent of the choice of 
the spectrum of turbulence ^(g), the emission intensity 
increases with frequency in the range of a; < ^LOp , and 
decreases when w > jcop , while the maximum is reached 
at cli = ^iOp. Then, instead of equation (^6|), we have 



P. = P 



A 



l^^l 



(38) 



The argument of this function has a minimum at a; = 
7Wp , and consequently the function achieves its maxi- 
mum at this frequency. However, we should note that in 
the case of convergence of the integral in equation (|3^), 
this maximum would be practically invisible. To demon- 
strate the behavior of P^^ at small frequencies, in figure || 
we show calculations for three different turbulence spec- 
tra ^ presented in the following specific form 



*(?) = 



Aai 



q2+a,^l^X^q2y-a,/2 



(39) 



Here, according to equation (py), the normalization con- 
stant 

"^~ r((l-ai)/2)' ^^' 

where T{z) is the gamma-function. The results of calcu- 
lations in figure correspond to three different values of 
Q!]^ = —1,0 and 1/2. It can be seen that while for 



Qf]^ = or 1/2 the integral in equation (^2|) diverges and 
the maximum of P^ is clearly seen at 70;^ , for the value 
of tti = —1 the emission intensity is characterized by a 
broad plateau without any distinct maximum. 

To explore the emission spectra in the frequency range 
jujp <C w <C Wj and w ^ Wp, let's assume that the 
turbulence spectrum has a broken power-law form: 



*(?) = 



A= 



/ \2+Qi 



AMf 






(41) 



where qi and 92 are constants of the order of 1/A , and 
the factor A"^ is introduced for the reason of dimension 
consistency. The condition for the convergence of the 
integral in equation ( p8|) on the lower and upper limits 
implies 

ai < 1 , a2 > 1 . (42) 

Depending on the value of ai there are two different 
cases related to the convergence of the integral in equa- 
tion ( p2| ) . If the integral is converging at the lower limit 
(i.e., ai < 0), we have the case discussed above. Let's 
consider now the range oi < ai < 1. Then, for the 
frequency interval ^iOp <^ lo <^ loj we have 



P, 



where 



e^X^qfiB^ 



Ci 



/2c7V 
-I- 3ai 4 



a 



■al 



(3 + ai)(2-Kai)- 
At lower frequencies, ujp <^ lj <^ l^p j 



— - (43) 



(44) 



Pu,= 



i^\\l{B^) 
27r2 m^c^ 



f2cLoqiY' 



- 1 



-• (45) 



9i 

allow a smooth 



We note that these equations ( [43| , | 

passage to the limit a^^ — > 0. 

In the range of large frequencies, a; S> Wj , the radiation 
spectrum behaves as a powe r-law uj~"^ , i.e. mimics the 
turbulence spectrum (TF87), 
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Fig. 2. — Spectral power calculated for the turbulence spec- 
trum ^ given in the form of equation (M) for three different 
values of power- law index a^ . 



where 



C-2^ — 



a2 l + a2 (2 -I- a2)^ 



^2 



(47) 



The energy lose rate of a charged particle due to radia- 
tion in t he magnetic field is given by the classical formula 
(see e.g. Landau & Lifshitz 1975| ) 



/ = 



2eV 



(/3 X Bf = 






B"^ sin^ X ■ 



(48) 



where x is angle between the particle velocity and direc- 
tion of the magnetic field. By averaging J, first over the 
directions, then over the strength of the magnetic field, 
and taking into account that (sin x) — § j ox\e finds 



/ = 



4eV 



IB' 



(49) 



(in the numerator, 13"^ — \ is substituted). By definition, 
the same result can be obtained by direct integration of 
equation (^ or equation ( p3[ ) over the emitted photon 
frequencies: J^Pujdoj = I. Nevertheless, it is worth to 
perform we such computations; they can serve as a good 
test for the consistency of the results. 
From equation (33) we find 



Pujduj = 



(B' 



GTT^m'^c 



d^ 



u{^) / duj 



2cY 
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After the substitution of the new variable u ~ 
(2c7^)/^, the integration over w leads to equation 
and then we obtain 



P,.,dLJ 



2eV(B' 



Zm 



2„3 






(50) 



(51) 



The remaining integral is equal to 2/3, so the direct inte- 
gration of the emission spectrum leads to equation (|49|). 
This interesting result, when the integration of the ap- 
proximate equation (^) gives precise expression for en- 
ergy losses, has a quite natural explanation; while equa- 
tion ( |33| ) contains the first (quadratic) term in the ex- 



pansion of the spectrum over the magnetic field strength 
B , the precise expression for the energy losses given by 
equation (H9) is proportional to the second power of B . 
Also we note that the energy losses are independent of 
the spectrum of turbulence ^ . 

In a similar way one can find the angular distribution of 
the emission after integration over the frequencies; there- 
fore we simply write down the final result: 



dir, 



4eV(B^) 1 + 7' 



,4/14 



37rm2c3 (1-^726*2): 



■ dn. 



(52) 



Here, din is the energy emitted into the solid angle dQ 
per time unit. This angular distribution also does not 
depend on the turbulence spectrum ^ . 

4. LARGE SCALE TURBULENCE 

In the case of a large scale magnetic field turbulence, 
A > Rl , the conditions imposed by equation (|^ are vio- 
lated, therefore the results of the previous section are not 
anymore valid. On the other hand, the radiation spec- 
trum formed in the regime A 3> Rl can be derived ana- 
lytically. In this case the particle deflection angle exceeds 
1/7, and the radiation spectrum, P^j , is determined by 
the instant curvature of trajectory (or the instant value 
of the magnetic field) . Thus the result should be similar 
to the spectrum of syn chrotron radiation in the homo 
geneous nragn e tic field (pchwingci 



[194'4: pinzburg fc Sy 
It the charged 



rovatskii 196£ ; Landau fc Lifshitz| 

particle moves perpendicularly to the magnetic field, the 

emission spectrum is determined as 



PUt) 



V3e2 
2ttRl 



Fix); 



where 



F{x) =x I K^/y,{u)du 



(53) 



(54) 



Here K^i^{u) is the modified Bessel function, x = uj/lOc 
and Wc is determined by equation (O). If the charged 
particle moves at an angle x to the magnetic field, in 
equation (p^) B should be substituted by th e perpen- 
dicular component of the field, B±^ = Bamx (Ginzburg 
fc Syrovatskii||1969p 

It the magnetic field is turbulent, then the spectrum 
Pi^{t) should be averaged over directions of the field, i.e., 
integrated over the pitch angle \ . This results in the 



following expression (Crusius & Schlickeisei 1986) 



P^{t) = 



V3e2 
2ttRl 



Gix); 



(55) 



where 

G{x) = — [W.^iix] W,^.{x) - W.^i{x)W_i^i{x)) . 

(56) 
Here Wfj,,a{x) is the Whittaker function. 

The function G{x) can be presented in a more conve- 
nient form: 

G{x) = ^ [(^ + 3a;^) (ni/sf + X K2/3 (2ki/3 - 3x^2/3)] , 

(57) 



Kl/3 



via familiar Bessel functions 
K^jzixjl) (lAharonian et alj pO10|) 
asymptotic befiavior both at 



G{x) has a simple 
ow and high frequencies: 



xRf^{x) 



G{x) 



^(r(i/3))V/^ 



2L p-^ 
2 



X< 1, 
X> 1. 



(58) 



Although differences between the spectra of syn- 
chrotron radiation in homogeneous and (large scale) 
chaotic fields, i.e. between functions F(x) and G{x) , are 
not dramatic , vet thev not too sma ll to be neglected in 
calculations ( Aharonian et al. 2010 ). In particular, these 
functions achieve their maximums, max(i^) — 0.918 
and max(G) = 0.713, at different points, x — 0.286 
and X = 0.229, respectively. Obviously, similar differ- 
ences we should expect for the spectral energy distri- 
butions described by the functions x F{x) and x G{x) . 
Namely, max(xF) = 0.693 is achieved at a; = 1.33, and 
max(a;G) = 0.444 is achieved aX x — 1.15. 

Finally, we note that function G{x) can be approxi- 
mated by a simple analytical expression. 



G{x) 



1. 



ci/3 1 + 2.21x2/3 + 0.347 a;'*/^ 



Vl + 3.4x2/3 1 + 1.353x2/3 + 0.217x4/3 
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whic h prov ides better than 0.2 % accuracy (Aharonian 
et al. 201C). Thus, this approximation can be safely used, 
instead ot the precise equation (|57|), in detailed calcula- 
tions of radiation in environments with large scale tur- 
bulent magnetic field. 

When deriving Eq.(Ba) we assumed that the magnetic 
field B is oriented chaotically, but its absolute value, 
|B| , is fixed. However, in a turbulent medium the spatial 
variation of the field strength could be quite significant; 
therefore we have to average the r esults also oyer th e 
absolute value of the field (see, e.g., Bykov et al. 2012). 
Let's introduce the distribution function 



w{B)dB = h„{b)dB/Ba, 



(60) 



where Bq = J (B^) , b = B/Bq . By definition, w{B) dB 

is the probability of the strength of the magnetic field 
being in the interval {B,B + dB) . We will consider three 
different distributions with the function /i„(6) presented 
in the following forms: 



hi{b) = 
h2{b) 



ho{b) ^ S{b - 1) , 
3^/6 



,^2g-3bV2, 



32 &2 



7r(l + 62)4 



For all distributions 



hn{b)db^ b'hn{b)db=l. 





(61) 
(62) 

(63) 
(64) 



For the variance of these distributions, D = {b"^) — (6)2 = 
1 — (6)2 , we have 
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Fig. 3. — Spectral energy distribution of synchrotron radia- 
tion corresponding to distributions of the magnetic field ho (solid 
curve), h\ (dashed curve), h2 (dotted curve). 

The energy lose rate of particles given by Eq.(^) de- 
pends only on {B ) , therefore it is convenient to compare 
the averaged spectra for the same value of (B^) : 



CX) 

(P^) ^ J Pu w{B) dB = IRhix)/u;o , 



(66) 



Do = 0, Di= 0.15, £>2 = 0.28. 



(65) 



where ujq = 3ei?o7^/(2TOc) , x = cj/cjq, /q Rh{x)dx = 
1. 

In Fig ra we show the SED of synchrotron radiation, 
xRh{x) , calculated for magnetic field distributions given 
by Eqs.(|l|),(|2|) and (|3|). 

One can see that the spectrum of synchrotron radia- 
tion is somewhat shifted, depending on the distribution 
of the magnetic field, towards higher energies, and, more 
importantly, it is significantly broadened compared to 
the radiation spectrum relevant to the (5 -functional dis- 
tribution of magnetic field hg . One can show (see Ap- 
pendix |A|) that if w{B) is characterized by a power-lay 
asymptotic dependence: w{B) ex B~'^ for B — >■ c», then 
the spectrum of synchrotron radiation also has a power- 
law asymptotic, namely (P) oc a;~'^+2 for u ^ ^q. Note 
that the power-law index should exceed cr > 3 in order 
to provide convergence of {B ) . 

It is important to emphasize the broadening and the 
shift of the spectrum of synchrotron radiation in a large- 
scale turbulent field has a quite different origin and 
should not be confused with the effects related to the 
jitter regime of radiation in a small-scale turbulent field. 
In this regard we should note that a broadening of syn- 
chrotron radiatio n has been "observe d " in the numeri- 
cal simulations of [Teraki fc Takahara] ( ^011 ). However, 
most likely the authors misinterpreted the obtained spec- 
tral feature and refereed it to an intermediate regime be- 
tween the synchrotron and jitter regimes of radiation. 
However, we believe that this component of radiation re- 
vealed in their simulations, has a standard synchrotron 
origin, but simply broadened because of the distribution 
of the strength of the magnetic field. 



5. ENERGY SPECTRA OF RADIATION IN THE 
JITTER AND SYNCHROTRON REGIMES 

In this section we compare spectra of synchrotron and 
jitter radiation, produced in two different large- and 
smaU- scale turbulent magnetic fields but with the same 
value (B ) , thus the total radiation power given by equa- 
tion (^) is the same for radiation in both regimes. For 
comparison, it is convenient to introduce the normalized 
the emission intensity: 



xR(x) 



R{x) dx — Pcj duj/I , 



: Uj/Ujc 



(67) 



Obviously, the following condition is held: /^ R{x) dx — 
1. 

For calculations we have to select a assume a spectrum 
of turbulence and a distribution of the magnetic field 
strengths. For the sake of simplicity, below we consider 
the case of chaotic synchrotron emission, i.e., we adopt a 
field distribution corresponding to equation (pl[). Then 
the function R depends only on the magnetic field. 



Ri-)-'-^Gix) 



(68) 



We note however that i? as a function of argument x 
does not depend on (B^) . In this section, for calculations 
of jitter radiation we consider the turbulence spectrum 
as an one-parameter family of functions: 



vl'(g) 



X^Ar 



(l + A2g2)l+a/2 



(69) 



The normalization constant, Aa , is obtained from equa- 
tion (H): 

_ 8.3/2 r(i + ,/2) 
"" r((.-i)/2) ■ ('") 



The spectrum presented in the form of equation (|69|) is 
characterized by a power-law dependence for q ^ A~^ . 
Although the spectra of turbulence, which can be gen- 
erated in astrophysical environments, remains an open 
question, usually it is approximated as a power-law. This 
assumption is justified by a few fundamental consider- 
ations. In particular, the power-law spectra of turbu- 
lence with spectr al indices of 5/3 an d 3/2 appear in the 
hydrodyna mical ( Kolmogorov 19411) and mag netohydro- 
dynamical (Iroshnikov 1963; Kraichnai^ 1965) turbulent 
media. 



Note that the asymptotic form of equation (69) is con- 
sistent with equation ( pl| ) for a^^ = —2, a2 = a. In 
figure y we show the normalized spectral energy distri- 
butions of the synchrotron and jitter radiation, xR{x) 
produced by particles of fixed energy jmc^ . The spec- 
tra are plotted as a function of a; = e/ec = huj/hujc , 
for three different indices characterizing the turbulence, 
a = 2,5/3,3/2. 

For a rather broad range of variation of the index a , 
from 3/2 to 3, the presentation of the turbulence spec- 
trum in the form of equation ( p9|) allows simple analytical 
approximations for the radiation power 



Pujduj = I f(xj) duj/iOj , 
where Xj — u) /loj , and 

f{xj) = Ca (1 + 0.22 Xj + 0.43 X 



2\-a/2 



(71) 



(72) 
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Fig. 4. — SED {xR{x) , see equation (pn)) of radiation of mo- 
noenergetic particles in turbulent magnetic field in the syn- 
chrotron and jitter regimes: synchrotron (black line) and jit- 
ter (red, green and blue solid lines). The spectrum of tur- 
bulence was taken in the form of equation (69). The ra- 
tio of the correlation length to nonrelativistic Larmor radius 
{Rl = mc? /eB) was adopted to be X/Rl ~ 10~^ . Red, 
green and blue lines correspond to the indices a — 2, 5/3 
and 3/2, respectively. 

The coefficient C^ is determined from the normaliza- 
tion L f(xj)dxj — 1. The comparison with the exact 
numerical calculations shows that the precision of this 
approximation is better than 7%. 

Figure demonstrates the basic spectral features of the 
jitter radiation. The SED peaks at energy which com- 
pared to the maximum of the synchrotron radiation at 
e — 1.155ec is shifted by the factor of 2/3Rl/X. Below 
the maximum xR{x) ex x, i.e. the SED increases with 
energy slower than the SED of the synchrotron radia- 
tion, xR{x) ex cc"*' 3 . Moreover, while the synchrotron 
spectrum has a quite sharp (exponential) cutoff beyond 
a; ^ 1 , the SED of jitter radiation after the break at 
x ^ Rl/X continues as a power-law, xR{x) oc a;^~" up 
to j:~ {RL/Xf. 

In astrophysical environments, acceleration of particles 
typically leads to broad energy distributions. Below we 
compare the synchrotron and jitter radiations for differ- 
ent distributions of accelerated particles A^(7) : 



oc 

/ 



F(w) = / P^ iV(7) d-1 . 



(73) 



Here we assume that energy distribution of all particles 
occupies certain energy interval (7min,7maa:)- Outside 
this interval, the function N is null Q For the jitter 
radiation, using equation ( |33|) and introducing a new di- 
mensionless function '^i{Xq) — '^{q)/X^ , as well as sub- 
stituting the integration variable 7 by ry = Aa;^/(2c7^), 

^ We would like to iad icate to the non-physical lower limit in the 
integral in equation (|73|). However, this convenient for integration 
representation is correct as long as the function N is taken zero 
outside the physically meaningful region 
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we obtain 



Piu) 






d^uiO 



^W$*'C"K\/^ 



^ • (74) i 



Let's assume now that the relativistic charged particles 
have a power-law distribution, N^j) = iVo7~'' . It can be 
shown that for the range of the power- law index, 1 < /i < 
2a-|-l , the main contribution to equation (|73| ) is provided 
by particles of energy 7 ~ (Aw/c)^/^. Therefore, for the 
energy interval u) ^ c/A , equation ( [74|) can be integrated 
over drj in the limits from to 00 : 



P(^) 



*A(B') /2c\ 



(A'-1)/2 



12 7r2m2c4 VAwy 



X /dfu(Or^^+'^/' /d?7*i(?7)77(''+^)/2_ (75) 



The power-law dependence of the spectra (-P(w) ~ 
[^;-(a'-i)/2 ~) jg explained by the same reason as in the case 
of the synchrotron radiation: cj and 7 enter into P^^ in 
a combined form w/'y' ^ (for discussion of the case of syn- 
chrotron radiation see Rybicki & Lightman 1979). Thus, 
for a power-law particle distribution the synchrotron and 
jitter mechanisms lead to the same type of energy spec- 
tra, therefore the ratio of the emission intensities due to 
these two processes, 



-* synchr vi-^ j 



C{fi,a) 



t) 



(A'-3)/2 



(76) 



does not depend on photon energy w . Interestingly, the 
index of /i = 3 appears to be special, in the sense that 
independently of the turbulence spectrum, the ratio r = 
1 . This, in particular, can be seen in figure at low 
energies. Note that although the energy losses due to the 
synchrotron and jitter mechanisms in the large and small 
turbulent fields are equal (for the same mean magnetic 
field), formally for fi > 3 larger energy is radiated out 
due to the jitter mechanism ( r > 1 ) , and vise versa, r < 
1 for /z < 3. This apparent inconsistency is related to 
the assumption of pure power-law particle distribution. 
However, for a realistic distribution of particles with a 
high energy cutoff, the spectral shape of the synchrotron 
and jitter radiations differ significantly. In particular, for 
power-law distributions with an exponential cutoff, given 
in a rather general form 



iV- 



7 ^ exp I 



il/lc 



\t^\ 



(77) 



in the high energy limit, the shapes of the synchrotron 
and jitter radiations spectra differ significantly (see fig- 
ure q). While the synchrotron co mponent beyond the 



ma ximum dec r eases exponentially (Lefa et allpOlS , also 



see Fritz 1989| ; ^irakashvili fc Aharonian 2007| 
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Fig. 5. — SED of synchrotron (dashed lines) and jitter 
(solid lines) radiation calculated for electron distribution: 
7~'' exp [— (7/7cut) ] • Calculations are performed for two val- 
ues of /3: /3 = 1 and /3 = 4 (the number labels indicate 
the used values for different lines ). The cutoff energy is set 
7cut = 10* , and the computations are performed for B = 1 G. 
The jitter radiation is computed for turbulence spectrum in 
the form of equation (BOl) with a — 5/3 and the ratio of the 
field correlation length to the nonrelativistic Larmor radius 
{Rl = mc^/eB) is \/Rl =0.1. 



with 



Wcut 



3eB 

2mc 



"fcut I 



(79) 



the jitter emission spectrum beyond the break around 
Wcut(^L/A) , has a long power-law tail, P(w) oc uj^" , in- 
dependently of the shape of the particle distribution in 
the cutoff region (i.e. the value of /3). In this regard, 
this is a unique feature of the jitter radiation which pro- 
vides direct and mo del- independent information about 
the spectrum of turbulence. As long as the condition of 
small-scale turbulence is satisfied ( A < Rl ) , we should 
expect radiation with characteristic a broken power-law 
type spectrum. While the photon index at low energies 
is directly related to the spectral index of relativistic par- 
ticles, r — (^-1- 1)/2 , or in the case of a low energy cutoff 
or very hard particle spectrum below the cutoff energy 
(e.g. in the case of Maxwellian type distribution - see 
figure |6|) , r = 1 , the spectrum after the break depend 
only on the spectrum of turbulence. In an environment 
with large scale turbulence, the picture is just opposite. 
The radiation proceeds in the synchrotron regime and 
therefore is not sensitive to the details of the turbulence. 
On the other hand, the synchrotron radiation carries in- 
formation about the overall spectrum of parent particles, 
including the most important (from the point of view of 
the acceleration theory) cutoff region. 

In Fig. we show a comparison of the synchrotron and 
jitter radiation for the case of Maxwell distribution of 
emitting particles. 

6. ANISOTROPIC TURBULENCE 

If the distribution of the charged particles is isotropic, 
the analytical solutions derived in the previous sections 
can be generalized to the case of the correlation tensor 
with an arbitrary angular structure. Indeed, similarly to 
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Fig. 6. — The same as in figure H, but for Maxweliian- 
type distribution of ciiarged particies: 7^ exp(— 7/7cut) with 
7cut = 10^ . 



equation (|27|), the radiation power can be expressed as 



Pr^Ut) 



72(1+7^04) 



{Spa- -VpVa) 



n'^m^c^ (1+72612)4 ^ p" p 



Kpa{q,x)2TT5{Cb + cqu 



. (Pa dx , , 
(iTTJ-^ ivr 



To obtain the radiation spectrum, this equation should 
be integrated over the photon emitting angles and av- 
eraged over directions of velocities of emitting particles, 
V = /3/|/3|- In order to simplify calculations, let's intro- 
duce the following intermediary tensor 



T. 



pa — 2 I ' P*^ 



'^p^a) 



jHi + j^0^) 



" "' (1-^726*2)4 

X(5(a) + cqv — yc) (Kl (Kla , 



(81) 



where dH. and dfla are the solid angles related to the 
directions of momenta of the emitted photon and the 
emitting particle, respectively. Note that equation (pm 
contains all the "directional" terms. Thus, the radiation 
spectral power can be expressed as 



P^it) 



TpcrKpcr{q,>c) 



d^q dx 



(27r)3 27r 



(82) 



The correlation tensor Kpa^ and its Fourier transforma- 
tion Kpa- are defined in Sec. 3. 

According to equation (p^), the tensor Tpa- is a sym- 
metric tensor of the secondorder; it depends only on the 
vector q , therefore, Tp^ has the following structure: 



Tpa — FiSpa + F2qpqa , 



(83) 



where -Fi and F2 are scalar functions. The convolu- 
tion of tensors Tpa- and Kp^ , taking into account the 
transversality condition of equation (UW , gives 



-^ per ^ pa — ^ 1^^ pp : 



(84) 

This expression determine the integrand in equation (|82|). 
To obtain the scalar function Fi the following relations 
can be used: 

Tpp = 3^1 + q^F2 , qpqaTp^ = <z'^l + q^F^ , (85) 



which give 



-^1=2 i^PP ~ IpqaTpa/q'^ 



(86) 



Using equation (pll), one finds 



1 



Fi = - 1 



(gz^ 
g2 

X (1+^2^2)4 ^(^ + ^1^ - ^) '^^ ^"" ' (87) 

which after the integration can be presented in the form 
(see Appendix 0) 



F.^-u{txm. 



Here U is determined by equation ([B5[). 

The trace of the correlation tensor can be represented 
as 

kpp{q,x)^{B^)^{q,x), (89) 

where ^{q, x) satisfies the normalization condition: 



(90) 



d^q dx 

Thus, the radiation power can be represented as 
„ , , 6^(52) n d^q dx 

^-(^^-^y^^(^'"/^^)*^^'")(2^2^- 

(91) 
If the correlation function Kpc{r,t) does not depend 
on time, i.e. "^{q^x) = 2'ii5{x)'^{q) , then the integra- 
tion over X is trivial leading to 



^"(^) = 2£S^/"(^)*(^)®(^-i: 



q 



(92) 



In the derivation of this equation we took into account 
that 

[/(e,0) = u(C)e(e-l). (93) 

Here u(^) is defined by equation (pi), and Q{x) is the 
Heaviside step function (i.e., Q{x) = 1 if a; > and 
Q{x) =0 if a;<0). 

Obviously, in the case of isotropic turbulence, the gen- 
eral equation (£2) should coincide with equation (pi]). 
Moreover, equation (31) can describe even the case of 
anisotropic turbulence, if one substitute the function 
^((7) by the spectrum of turbulence averaged over di- 
rections of the vector q, vE'(q) = (^(q)). This implies 
that the averaged radiation power does not depend on 
the structure of the correlation tensor. In particular, 
the monotonic decrease of the intensity given by equa- 
tion ( pq ) , also is observed in the case of anisotropic tur- 
bulence. 

Although this conclusion is derived under the assump- 
tion of isotropic distribution of emitting particles, in fact 
the obtained result is valid also for the case of anisotropic 
particle distribution when the change of particle density 
is small for typical "angular shifts" of the value of I/7. 
Thus, as long as the emission is considered in the ultra- 
relativistic regime, this assumption can be broken only 
in the case of highly collimated particle beams (which is 
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likely inconsistent with the underlining assumption of the 
turbulent magnetic field). Also we note that if we deal 
with strongly anisotropic distribution of particles, the ra- 
diation does depend on the structure of the correlation 
tensor, therefore it is important to define it correctly (see 
Section 0). 

7. COMPARISON WITH PREVIOUS RESULTS 

In recent years, a large number of studies have 
been devoted to calculations of radiation (the magnetic 
bremsstrahlung) generated by charged particles in small- 
scale turbulent magnetic fields. However, to our knowl- 
edge, the general expression for the radiation spectrum 
described by equation (pl|), is derived for the first time 
in this paper. Also, in the previous studies a few ad- 
ditional conditions have been assumed, which however 
appear redundant, and actually not needed at all in the 
framework of our approach. This redundancy not only 
superficially constraints the applicability of the obtained 
results, but also introduces some confusion in the anal- 
ysis and comparison of different radiation regimes. Fi- 
nally, some solutions and related conclusions derived in 
this paper do not coincide with the results of previous 
studies. Therefore, we present below a short overview of 
a few important papers on the topic, compare their main 
results with our study, and outline the key differences 
between the approaches which might cause, in our view, 
these discrepancies. 

There are two basic theoretical approaches to study ra- 
diation in random magn etic fie lds. The first one is based 
on the seminal paper by TF87, where a kinetic equation 
has been derived for the probability of different particle 
trajectories in a chaotic magnetic field (see Eq.(12)TF87), 
and an approx imate solut i on ha s been found to this equa- 

2005| , for a simplified descrip- 

However, the introduced 

simplifirstinns signifirsntljf limit the applicability of this 



tion (see also Fleishman 
tion of the approach of 
ifiratinns signifi 



appr oach and do not allow a self-consistent testament of 
the problem. More specifically, we discuss these issues 
bellow. 

T he second approach is based on the pert urbation the- 
ory (|MO0| ; |Medvedev| |2006| ; |Fleishman|[2006| ). In all these 
papers, the authors start from an expression for the emis- 
sion produced by a particle defiected by a small angle in 
a m agnetic field localized in a compact region of space 
(see Landau fc Lifshitz 1975 , § 77). This expression can 
be written as 



du! 



a^ 






W'72 



2uj'^j 



'2-,4 



duj\ 



(94) 

where a^i = J_ a{t) e*" * dt is the Fourier component 
of acceleration. However, if the magnetic field occupies a 
large volume, then even in the case of chaotic magnetic 
field, the particle deflection will be (unavoidably) large 
(because of multiple scatterings over emission correlation 
length). Therefore, the solution based on this expression 
has a rather limited applicability compared to the prac- 
tical realizations in the chaotic magnetic field. 

Th e approached employed in our otudy alao io baood on 



the perturbation thenr^^; but it is ^fa1id when the partir 



tion length, or, equivalently, if i?L ^ A . 

Note that th e later approximation was also implicitly 
used in TF87 (see Eq.(ll)TF87) when deriving the ki- 
netic equation. So even in the case of precise solution 
of this equation, the results cannot be expanded beyond 
the parameter region described by the perturbation the- 
ory approach presented in our paper. Moreover, since 
the derived kinetic equation appeared to be too complex 
to be treated analytically, a few further simplifications 
have been introduced to obtain an analytical solution. 
In particular, the original Eq.(12)TF87 was replaced by 
Eq.(17)TF87j which indeed could be equivalent to the 
original one if in the rhs of this equation they would 
use the function q{uj,9) determined by Eq.(15)TF87- In 
Eq.(17)TF87 "^ enters as a parameter, thus for solution 
of this equation it can be taken as a constant, and the 
function q treated as a function of one variable q(9) . 
However, since in this case the equation does not have a 
solution, the authors replace the function 5(0;, 6) by an 
empirical function q{uj) . This simplification allows an 
analytical solution, but since it concerns the term with 
the highest derivative in the equation, the uncertainties 
imposed by this substitution cannot be evaluated and 
correspondingly, the limits of applicability remain highly 
unknown. Note that the empirical function q{uj) itself 
determines the radiation spectrum in the case of absence 
of the regular component of the magn etic fie ld. How- 
ever, within the framework of theory of TF87| this func- 
tion, strictly speaking, is not derived. Instead, based on 
arguments of the asymptotic behavior, they have pro- 
posed the following form q{uj) — q{uj,9 — 0^,), where 
^* = (a — l)/7^ and the value of parameter a was de- 
termined "from the requirement that at high frequencies, 
where the perturbation expansion (the method of equiva- 
lent photons, see Appendix) is valid, the present method 
yields the same result as the pertu rbation expansion" 



Let's consider now the results of Fleishman & Bieten- 



holz (2007), where the approach of [TEST] has been ap- 
plied to the case of the random magnetic field without a 
regular component. To make the comparison transpar- 
ent and less bulky, discuss the results for the fix value 
of the index of the turbulence spectrum a = 2 , and ig- 
nore the impact of the surrounding medium (i.e., assume 
uJr, — 0). For this specific case, the spectrum obtained in 



Fleishman & Bictenholz (2007), can be expressed as: 






Sttc 



Here 



q{uj) 



^It^ol"^ 



{auj/2Y + (c^o7')2 



00 

$(s) = 245^ I dt e-^"* sm{2st) (cotht- - 



(95) 
(96) 
(97) 



where 



1 



872 yqiuj) 



1/2 



c c 

t^st = t;- , Wn = - . (98) 



Rl 



A 



$(5) has the fol lowing asymptotic limits (Fleishman 



_Bietenhoy gOOTD 



deflection is small on the typical magnetic field correla- 



$(s)wl, ifs>l; $(s)sa6s, if s < 1 . 



(99) 
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At a; ^ Uj — ujfff^ the s parameter is large, s ~ 
Uq/uOsi — Rl/^ >• 1- Thus one can use the asymptotic 



Hmit for $(s) = 1 
in a simple form 



Then equation (B5[) can be expressed 



P.., 



3eV 
Sttc 



9(^): 



(100) 



Apparently, the function q{uj) determines the shape of 
the radiation spectru m. However, this function has not 

( |2007| ) 
thi 



be derived either by Fleishman & Bietenholz 
by Fleishman ( 2005 ). We can only guess that t 



or 
le authors 



have used the simplified form of Eq.(39)TF87 (after re- 
moval of the bulky complex term from that equation) . 
' Fleishman & BlcLeuholii' (2007) performed numerical 



calculation of the radiation spectra also for the case 
of strong random magnetic field, i.e. in the regime of 
A > Rl ■ In the asymptotic limit of A >• i?L , the spec- 
trum can be obtained analytically; in this regime we deal 
with the standard synchrotron spectrum described by 
equation (55). On the other hand, equation (By) with the 
asymptoticTimit of ^(s) for s ^ 1 from equation I^M), 



differs significantly from equation (55|). In our view,~the 
reason for this discr epancy is that the basic kinetic equa- 
tion in the theory of TF87 is derived under assumption of 
A ^ Rl (see Eq.(9jTF87)- Thus, this approach cannot 
be applied for the regime A > i?L • 

The fact that equation (BSh is not applicable for the 
case of A ^ Rl can be also illustrated by computation of 
the total power emitted by a particle. Let's consider the 
ratio of the radiated and lost energies by the relativistic 
charged particle: 



P(Ld) did , 



(101) 



Here / and P^ are determined by equations ( p9| ) and 
(|95[), respectively. For a particle emitting in vacuum, the 
condition p = 1 should be satisfied. I n th e asymptotic 
case of A ^ Rl , one can use equation (|lO0| ) and demon- 
strate that for a = 2 we indeed have p = I . However, in 
the limit of A » Rl , 



p = 5.6 X p 



2/3 



where 




(102) 



(103) 



i.e. the condition p = 1 is violated. 

Thus, we can conclude that in the case when the non- 
cha .otic m agnetic field is nil, the approached developed 
by TF87 has a very limited applicability. Namely, one 
can derive the spectrum in the form of equation ( |100| ) 
with function q{Lj) constrained by asymptotic behavior 
only. 

Now let's compare our results with the studies based, 
like our pa per, on the perturbation theory. 

In M00| , a specific geometry of interaction has been 
postulated. Namely, it was assumed that the particle 
moves along axis x, and that magnetic field has only y 
component, By . Therefore, acceleration is parallel to z- 
axis: ait) = ^ By{vt,0,0) = ^B{t). Although the 
radiation power obtained in Sect. ^ was derived under 



the assumption of homogeneity of turbulence , and thus 
is not applicable to the case considered by MOO , it is 
straightforward to apply our approach to this case also. 
Namely, accepting the definition of the spectral power 
given by equation (0) , one can average over the magnetic 
field configurations m equation (^4|). This gives 



P^it) 



{B^ 



,,'2 



1 



27rm^7^c^ 



U} U} 



duj', (104) 



u/{2^^) 



/(27=) 

where 

oo 

iB^U= f{B(t + T/2)B{t~T/2))e'^'^dT. (105) 



is the Fourier component of the magnetic field correlation 
function; generally it may depend not only on w' , but 
also on t . Obviously, 



{B'it))=J{BX^ 





(106) 



Equation (104) describes the radiation power for the 
geometry postulated by for an arbitrary spectrum of tur- 
bulence. Then, a s pecifi c spectrum of turbulence has 



been considered in MOO for derivation of Eq.(17)Moo- 
The interesting feature of the latter is that in the limit 
of cj — > the spectrum P^^ ^ cj , and contains an 



abrupt cutoff, P^ 



at 



> 2a;,- 



However, we 



should note that these spectral features do not charac- 
terize the jitter radiation in general (as it can be seen 
from equation (104)), but simply are the consequence 
of the choice of a specific turb ulence spectrum and /or 
interaction geometry (see also Fleishman 2006| ). In- 
deed, if one adopts a different turbulence spectrum, e.g. 



(B2 



(o; 



/2 



oj: 



':) " , then for any positive value of 



the index a , the spectrum is a monotonically decreasing 
function of uj . Moreover, if the stochastic field has both 
y- and z- components, and the correlation function is az- 
imuthally symmetric in respect to x-axis, then ev en for 
the power-law spectrum of turbulence, adopted by MOC , 
the spectrum is not expected to be linear in the limit of 
small id . 

The treatment of radiation in a chaotic magnetic field 
always leads to the appearance of the correlation ten- 
Kf^a (see eq uation ( [19|) of this paper and Eq.(12) of 
However, often the structure of this 



sor. 



Fleishman 200' 



eq u 

3- 



tensor, Kpa-{q, x) , is wrongly postulated. If we consider 
a homogeneous environment without preferred direc- 
tions, then we deal with only two second-order tensors: 
5pa and qpqa ■ Therefore, the correlation tensor should 



= Ci 5pcr + ci qpqa , 
The transversality 



have the following structure: Kpa 
where ci_2 are two scalar functions 

condition implies ci +C2(7 

portional to [bpa — qpqcr/q^), as used in equation (|2^). 
However, in the case of existence of a distinct direction, 
s , e.g. normal to the shock front, the tensor structure 



2 — , thus Kpa has to be pro- 
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becomes more complex: 

Kpa = Ci 5p^ + C2 qpQa + Cs SpS^+Ci qpS^ + C^ Spqo- ■ (107) 

The transversality condition imposes three constraints on 



the scalar functions 



thus the correlation tensor K 



must have the following structure: 



per 



K, 



+ *2 Sp - qp 






qpqa 



'pa 



qa 



(sq) 



(108) 



where functions \I'i^2 depend on three arguments: |<7|, 
(sq) and >c. In a gyrotropic medium, the correlation 
function may contain an additional term: ep^^q^"^^ ^ 
where ep^r is the Levi-Civita-tensor and ^3 is a com- 
plex function as it fo llows from the general t heory of 
fluctuations (see e.g. Landau fc Lifshitq 1980| , § 122). 
Note, however, that this term does not contribute to the 
emission power, since in equation (g^ the tensor Kp^ is 
convolved with a symmetric tensor. 

For the additional assumption that the magnetic field 
is perpendicular to the direction s , the correlation func- 
tion should satisfy the equations Kpa-Sp = , Kp^Sa = . 
In this case the functions ^1 are linked ^2 via the rela- 
tion 



*i + (l-(sq)7g^)*2 = 0. 



(109) 



and the correlation function Kp„ is determined just by 
one scalar function. 

However, in some studies dealing with the anisotropic 
turbulence, different tensor structures have been pro- 
posed for the correlation function: Kp ^ oc (5nr, ^^nSg) 
- we can re fer, for example to Eq.(8) in jMcdvcdcv 2006 , 
Eq.(1 8) in |Fleishman| ^006[ Eg. (10) in [Mcdvcdev cTal 
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Onl , Eq.(ll) in [Reynolds fc Medvedev|p012| ). This cor 
relation function does not satisfy the transversality con- 
dition, i.e., the considered magnetic field is not diver- 
gence free, VB ^ . Apparently, this is a wrong result, 
therefore the results obtained under the assumption of 
Kpa (X (5 pa — SpSa-) should be revised. 

It is important to note that certain mathematical op- 
eration often used for computation of emission in chaotic 
magnetic fi eld lack matheniatical strictness ( also see dis- 
cussions in Fleishman 2006; Mcdvcdev 2005). In partic- 
ular, this concerns the involvement of the field harmon- 
ics, Bk , which implies that the Fourier transformation 
can be applied to the stochastic magnetic field. This 
assumption hardly can be justified or disproved from 
the first principles, however this assumption may lead to 
a rather controversial expression for the Fourier image 
of the correlation function. For example, the following 
structure has been obtained for the correlation tensor 
Kp„ (X B, 
man 



2006 



Brr fsce footnote 2 and equation 5 in Fleish- 



Medvedev 2006, respectively), which, ho' 



ever, cont radi cts the g eneral tensor structure g iven by 



equation (10S)(see also Landau fc Lifshitz 198C, § 122) 



Finally, we note that in the framework of our approach, 
no assumptions regarding the properties of the stochas- 
tic magnetic field are required. Instead, we assumed that 
the Fourier transformation can be applied to the mag- 
netic field correlation function, which is a significantly 
less demanding assumption. 



8. DISCUSSION AND SUMMARY 

The so-called jitter radiation mechanism represents a 
version of the magnetic bremsstrahlung of relativistic 
charged particles in a turbulent magnetic field. This 
regime of radiation can be realized with an efficiency as 
high as the "standard" synchrotron radiation but with 
quite distinct energy spectrum strongly shifted towards 
higher energies. This makes the jitter radiation an at- 
tractive gamma-ray production channel in highly turbu- 
lent astrophysical environments. 

In this paper we present a novel study on spectral prop- 
erties of the jitter radiation performed within the frame- 
work of perturbation theory in the regime of the so-called 
small-scale turbulence, when the coherent length of the 
field is significantly smaller than the nonrelativistic Lar- 
mor radius, A <C Rl = mc^ /eB , or 



A < 1.7 X 10^{m/me)iB/l G)" 



(110) 



Here B is the average magnetic field, and m is the mass 
of radiating charged particle. It is remarkable that the 
condition for realization of the jitter regime does not de- 
pend on the particle energy, but only on its mass. For 
exam ple, for electrons the condition imposed by equa- 
tion ( 11C| ) implies a turbulence scale as small as 100km 
in young supernova remnants, less than 10m in gamma- 
ray production regions of blazers, and 1cm in GRBs, as- 
suming typical values of the magnetic field in these ob- 
jects of about 100/iG, IG and IkG, respectively. For 
protons these conditions are relaxed by three orders of 
magnitude. However, the magnetic bremsstrahlung of 
protons is a much slower process compared to electrons. 
It becomes adequately effective only at very high ener- 
gies of protons and at the presence of large magnetic field 
which in its turn implies tighter conditions on the turbu- 
lence scale. Whether turbulen t fie lds can be generated 
at scales imposed by equation ( |110| ) is a non-trivial issue 
the discussion of which is beyond the scope of this paper. 
Here we focused merely on the study of radiation proper- 
ties and p erfor m calculations under the assumption that 
equation (110) is (by definition) fulfilled. We derived an 
expression for the spectral power of radiation presented 
in a general but rather compact form convenient for nu- 
merical calculations. 

The jitter radiation has a simple spectral form. Its 
SED for a monoenergetic particle distribution is shown 
in Fig.Q together with the SED of synchrotron radiation. 
Both SEDs have pronounced maximums separated from 
each other by a factor of i?L/A. Obviously when the jit- 
ter regime is realized, the maximum of its SED is shifted 
towards higher energies (the position of the peak in the 
synchrotron SED is at the energy « 1.115wc). Unless 
one introduces strong assumptions regarding the turbu- 
lence spectrum and/or its geometry, the low energy part 
of the spectrum has standard photon index F = 1 . It 
is hard but still softer than the spectrum of synchrotron 
radiation with F — 4/3 . The jitter and synchrotron 
spectra are very different beyond their respective max- 
imums. While the standard synchrotron spectrum cuts 
off quite sharply (exponentially) just after the maximum, 
the spectrum of the jitter radiation continues as a power- 
law until the energy '^ {Rl/X)^uJc with a photon index 
F = a -|- 1 , where a is the power-law index of the tur- 
bulence spectrum (in the framework of the perturbation 
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theory, the spectral shape of radiation above this hmit 
cannot be calculated). Remarkably, this part of the spec- 
trum is not sensitive to the details of the energy distri- 
bution of particles, but depends only on the position of 
the cutoff in the particles distribution. The latter deter- 
mines the start of the power-law tail which should be (by 
definition) quite long since Rl 'S> A. For example, if the 
ratio i^i/A exceeds 10, the power-law tail of the jitter 
radiation, which mimic the turbulence spectrum, would 
span over more than two energy decades after the max- 
imum. Bellow the maximum, the jitter radiation does 
depend on the particle distribution. In particular, if the 
relativistic particles have a power-law distribution with 
an index //, the spectrum of the jitter radiation is also 
power-law with photon index T = {^-\- 1)/2 , i.e. exactly 
the same as in the case of synchrotron radiation. 

In this paper we do not aim to discuss astrophysical 
implications of jitter radiation which deserve separate 
consideration. Instead, we would rather comment on 
some interesting features which make this mechanism 
quite unique amongst other high energy radiation pro- 
cesses. 

First of all, the slight dependence (or rather indepen- 
denc e) of the high energy spectral tail on the distribution 



of p£ rent relativistic particles, is quite unusual and ap- 
parently does not have an analog in astrophysics. More- 
over, the fact that the spectral shape of this tail simply 
mimics the spectrum of turbulence, opens a unique op- 
portunity for the straightforward probe of the spectrum 
of small-scale turbulence by measuring the characteristic 
high energy electromagnetic radiation and identifying it 
with the jitter mechanism. 

While in the case of injection of relativistic electrons 
into a highlv turbulent medium, where the condition of 
equation ( |llO| ) is satisfied, guaranties production of ra- 
diation in the jitter regime, the questions of its detection 
depends on the total energetics in relativistic particles. 
Given the usually (very) high efficiency of jitter radia- 
tion, and for typical parameters characterizing the non- 
thermal astrophysical sources of both galactic and ex- 
tragalactic origin, the production of detectable fluxes of 
jitter radiation in the X-ray and/or gamma-ray bands 
could be readily realized in standard acceleration and 
radiation scenarios. 

The identification of the origin or radiation is the sec- 
ond critical issue. Fortunately, the jitter radiation does 
provide us with distinct features which can be used to 



identify its nature. In particular, for a standard power- 
law injection distribution of electrons with a high en- 
ergy cutoff given, for example, in the form of equa- 
tion ( [771 ) with ^ = 2 , and assuming a Kolmogorov- 
type spectrum of turbulence, a = 5/3, we expect a 
gamma-ray spectrum which can be described as bro- 
ken power-law. The high energy part of the spectrum 
above the break is expected to have a photon index of 
r2=Q: + l~2.7, while the low energy part (below the 
break) Fi = (/i -I- l)/2 = 1.5 or Ti^ {fi + 2)/2 = 2 for 
the slow and fast cooling regimes, respectively. This cor- 
responds to the change of the spectral index by AF = 1.2 
or 0.7 depending on the cooling regime. Such a behavior 
differs significantly from the standard synchrotron cool- 
ing break with AF = 0.5. In the case of a low energy 
cutoff in the electron spectrum, which is often required 
to fit the data, e.g. the spectra of gamma-ray blazars, 
we should expect another break below which the photon 
index would be F = 1 Therefore, in the case of detection 
of a non-standard broken power-law spectra, especially 
when the high energy power-law tail has a photon index 
close to 2.5 and extends well beyond the break, the jitter 
mechanism can be be treated as a proces s responsible 



for the observed spectral features (see also Fleishman 



Bictenholz 2007) 



Despite all attractive properties of synchrotron radia- 
tion of ultrarelativistic electron, its spectrum usually ter- 
minates before reaching the gamma-ray domain. Even in 
the extreme accelerators it cannot extend beyond the so- 
called synchrotron limit ~ 100 MeV, unless being ad- 
ditionally Doppler boosted in sources with relativistic 
Doppler factors. This can be the case, for example, of 
the recently discovered fl ares of the Crab N ebula (see. 



Buehler et al 2012; striani et al. 2013, and refer- 
s therein) or t he multrT:ie V counterparts of gamma- 



ray bursts (Abdo et al 



2009). On the other hand, the 



jitter mechanism may ofler another possibility for exten- 
sion of the spectrum well beyond the energy synchrotron 
limit. We should note in this regard t hat i n the case of 
fulflUment of the condition in equation (110), the appear- 
ance of jitter radiation is not only unavoidable, but its 
spectrum could extend to high or even very high ener- 
gies. A more a principal issue in this regard is the chal- 
lenge of formation of turbulence on very small scales, e.g. 
through the Weibel type instabilities. 



APPENDIX 
HIGH ENERGY ASYMPTOTICS OF SYNCHROTRON RADIATION 

At high frequencies, the intensity of synchrotron radiation decreases exponentially, G{x) ex exp ( — 3^"^ ) = 

expl— ^^^j (see equation (pq)). Therefore at oj 3> wq the contribution to the radiation spectrum given by equa- 
tion (|6^) is dominated by regions characterized by large magnetic field: 6 S> 1 . Let us assumed that for 6 S> 1 , the 
magnetic held strength is distributed as power-law : wdB « Ab~'^db. 

It is convenient to start the calculations of the spectrum from equation (m) , where magnetic fleld strength is replaced 
by _B — > Bsinx, i.e. before integration over directions of the magnetic fleld. This gives 



P:.= 



V3e^ 



2ttR 



LO 



Kz/z{u)du, 



(Al) 



x/b sin X 



where x — co/ujq^ Rlo — mc^/eBo and b — B/Bq. To average over the strength of the magnetic field strength, 
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one should multiply equation ( |A1[ ) to the distribution of magnetic field, Ab^'^db, and then integrate over db. Let us 
introduce a new variable 6 — )■ ^ = x/b sin x , and change the order of integration over dS, and du . Then, after a rather 
simple analytical integration over d^ , we obtain: 

P^wib) db^A ^^ x-^+\^mxr-^ fu^-^K,/,{u) du . (A2) 

ZttHlq cr - 1 J 



The remaining integral over du can be expressed in terms of Gamma functions: 



u'^-'K,^,iu)du = 2--^r(^- + -j r(^- - -j . (A3) 



The integration over pitch-angles also leads to an expression containing Gamma functions: 

((sinx)"-^> = ^^/^ ^ P . (A4) 

\v a; / 2r(f + l) ^ ' 

Thus, in the limit of large frequency to ^ loq, the spectrum of synchrotron radiation averaged over the directions and 
strength of the magnetic field, is described by a power-law function: 

{P^)=A^C,x--+\ (A5) 

where 

V3 2-^ r(| + |)r(f-|)r(f + i) 



Note that if the magnetic filed is formally distributed as pure power-law, B: w{B) — AB " , equation (A5) gives 
precise solution for the radiation spectrum. In this case A — AB^^"^ , and therefore (Pui) appears to be independent 

°f^o- .. ... . rn 

To understand the condition for applicability of equation ( [A5| ) , let us estimate the correction terms to this equation 

for a specific distributions of the magnetic field. Let's assume, for example, the following distribution: 
In the limit b ^^ 1, the first two terms of the series are 



(1 + 62)1+'t/2 yi^a 5^+2 

Correspondingly, 



{P^) =A^{C„ x-^+2 - (1-1- a/2) a+2 x-^) . (A9) 



The ratio of these two terms can be expressed as 

(cr-l)(9cr2_25) 



18x2 2a;2 



(AlO) 



which implies that one can neglect the second term when x > cr'^/^ This can be treated as the condition for 
applicability of equation (A5). 



THE CASE OF ANISOTROPIC TURBULENCE 



Here we present some intermediate calculations required for derivation of equation ( |88| ) from equation (87). To 
compute the integrals, it is convenient to introduce the following new variables: 

,^ = 72^,2^ a; = cos??, (Bl) 

where 1? is the angle between vectors v and q . The integration over the azimuthal angle is trivial; it gives dfl dfl^ — 

In' 



^'^ ■ dxdC,. Then, the integration of equation (pW) results in 



' dC / dx(l + x^) y S{Cj + cqx - x) . (B2) 



2J^J' ' (1 + C)^ 
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For the upper limit ol integration over C we set oo , which is vahd only in the uhrarelativistic regime (see also the 
discussion after equation (|30|)). The argument of the (5-function nulls for x — Xq — {x — uj)/cq. The integral becomes 
zero if Xq lies beyond the integration interval, Xg > 1 or Xq < —1. For the value of Xq within the integration range, 
i.e., I^ol < 1, we obtain 

„ C2 



Ci 



(B3) 



where the lower and upper integration limits, Ci^2 , are determined by the conditions JXqI — 1 and C ^ 0. 

It is convenient to express these limits as Ci = max(0,^(K — 1) — 1) and C2 = S,{n + 1) — 1 (C2 should be positive), 
where 

K = — , ^ = . (B4) 

qc 



This allows derivation of equation ( ^81 ) via the analytical integration: 

u{i. k) = e(e(Ac + 1) - 1) \Ur e(i - i[K - 1)) + u^ e(e(«: - 1) - 1)] 

where 

2C^ - 2^ + 1 C - 1 2 \ 3 

^Tl (k + 1)2 + («; + l)3y ~ ^ 

and 



C/i = ^(^('* + l)-l)(e'(«-l)+4e 



2 = TsT^ ^3 (4 + S'^'^' + 3Ck' - ee'K* - 12Ck3 + (Se" + 4) k2 + 9^k) - -3 (^k + 1) In - 



t/2 = 



(B5) 
(C«; + l)ln(e(^ + l)), (B6) 

(B7) 



Equation (B5) implies that function C/(^, k) has non-zero values only if k > 7- — 1. The two terms in equation 
C/1^2, give non-zero contribution for 7- — 1<k<7- + 1 and k > 7 + 1, respectively. The continues function U has a 
break at k = 7 + 1 C^(^, k) . 
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